Large diamagnetic persistent currents 
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Abstract 

In multichannel rings, evanescent modes will always co-exist with propagating modes. The evanes- 
cent modes can carry a very large diamagnetic persistent current that can oscillate with energy and are 
very sensitive to impurity scattering. This provides a natural explanation for the large diamagnetic 
persistent currents observed in experiments. 
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Biittiker, Imry and Landauer first suggested the possibility of observing persistent current 
in normal metal or semiconducting rings threaded by an Aharonov-Bohm flux-. This current 
is an equilibrium property of the ring, given by the flux derivative of the total energy of the 
ring. Since then several experiments have been done that confirm the existence of persistent 
current in such rings-^^^ 1 ^ through magnetization measurements. However, the nature 
of these currents are quite different from what is expected theoretically^^. While earlier 
experiments^^^ 1 ^ had some ambiguity, recent experiments made on an ensemble of 10 5 rings 
have made very careful measurements of the sign (positive implies diamagnetic and negative 
implies paramagnetic) and periodicity of the persistent current^. If an ensemble of rings is 
taken, one can calculate an ensemble average over the number of electrons in different rings 
or over disorder, or over both 10 . In such cases one finds that the ensemble average has <po/2 ( 
or hc/2e) periodicity and the low field persistent current is paramagnetic in nature. One can 
also take a fixed chemical potential and average over disorder. Here again one can calculate to 
find paramagnetic persistent current with (po/2 periodicity 11 . Whereas, the experiment^ shows 
a persistent current, that has 4>o/2 periodicity but diamagnetic in nature at low fields and also 
of a large magnitude (10 to 100 times larger than that theoretically estimated in the above 
mentioned models of ensemble averaging). 

Any quantity that is very sensitive to disorder will average to zero. But the second harmonic 
do not obey this rule and gives nonzero value. This is essentially because the second harmonic 
consists of time reversed trajectories and disorder configuration does not change the observed 
quantity randomly. This is very robust and manifests in a variety of phenomena briefly described 
below. Weak localization in disordered metallic or semi-conducting samples occur because of 
this. Forward scattering probability beyond a certain length turns out to be negligibly small 
while the back scattering arising due to time reversed trajectories always interfere constructively, 
irrespective of disorder configuration. As a consequence the Aronov-Altshuler-Spivak weak 
localization correction to conductance has 0o/2 periodicity 12 . Also the response of a long 
cylinder to an applied magnetic field turns out to have </>o/2 periodicity™ 0o/2 periodicity of 
ensemble averaged persistent current is due to the same reason that the first harmonic averages 
to zero while the second does not. 
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The first attempt to explain the discrepancy in sign and magnitude is based on repulsive 
interactions between electron a ^ ^ . This did not turn out to be the correct mechanism because 
this yields a paramagnetic response at low fields whereas recent experiment has conclusively 
shown that the observed response is diamagnetic. For a recent analysis of the effect of disorder 
and interactions, we refer [15]. A more recent attempt to explain the experimental discrepancy 
is based on additional currents that may be generated in rings due to the rectification of a high 
frequency non-equilibrium noise — . This mechanism can give a diamagnetic current in absence 
of spin orbit coupling and a paramagnetic response in presence of spin orbit coupling. Recent 
experiment 9 - has also ruled out this explanation as paramagnetic response could not be observed 
in presence of strong spin orbit coupling. The origin of a high frequency non-equilibrium noise 
also seems to be unclear. 

All experiments have been done at finite temperatures. There are thermal effects wherein 
an electron can get energy from some collision and get excited to higher states. Such inelastic 
processes will not destroy the persistent currents. Persistent currents are actually observed 
in networks, where total length is much greater than the inelastic mean free path^. When a 
mechanism for excitation is present then evanescent modes can be excited. It has been shown 
that in a one dimensional (ID) ring, evanescent modes can carry a diamagnetic persistent 
current. It has a very small magnitude compared to persistent current in propagating modes and 
cannot exhibit 0o/2 periodicity as it is not sensitive to disorder. Rings used in the experiments 
have a finite thickness and are referred to as quasi-one dimensional (Q1D) rings. In this work 
we show that in Q1D, evanescent modes can carry very large diamagnetic persistent current 
that are comparable to that of propagating modes and are as sensitive to disorder as that 
due to the propagating modes. So this mechanism is a natural explanation for the observed 
diamagnetic persistent current. 

In this work we use a simple technique to excite evanescent modes. We consider the ring to 
be coupled to an infinite wire as schematically shown in Figure HJ This basically constitutes an 
open system and it is known that it can simulate the effects of inelastic collisions and thermal 
effects 19 . We shall see in our mathematical analysis how evanescent modes are excited in this 
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system very naturally. We consider two modes of propagation as the results can be generalized 
to any number of modes. There is a ^-potential impurity present in the ring at any arbitrary 
position X [Fig [T] . We apply Aharonov-Bohm flux </> through the ring, perpendicular to the 
plane of the paper. The Schrodinger equation for a Q1D wire in presence of a (^-potential at 
x = 0, y = yi is 

l-^(^ + o^) + v c(y) + iS(y-yi)t{x)M^y) = E*{x, y ) 

(i) 

Here V c is the confinement potential making up the quantum wires in Figure [H Solutions 
to Schrodinger equation is a ring geometry can be obtained by applying periodic boundary 
conditions to Eqn. [TJ The magnetic field just appears as a phase of *$?(x, y) that will be 
accounted for while applying boundary conditions. Away from the scattering regions Eqn. [1] 
can be separated as 

_ h 2 d 2 ^(x) Ti 2 k 2 
2m dx 2 2m 

and 

r h 2 d 2 



■iP(x) (2) 



2m dy 2 



+ V c (y)]Xn(y) = E n Xn{y) (3) 



Here we take V c to be a square well potential of width W that gives Xn(y) = s ^ n [^(y + \ )]■ 111 
the first mode, k\ = yj — ^ and in the second mode k 2 = J — ^ are the propagating 
wave-vectors, m is the electron mass, E is the electron energy and W is the width of the 
quantum wire. When electrons are incident along region / (in Fig [T]) in the first mode the 
scattering problem can be solved exactly. The solution to Eqn. [2] in region / becomes 

fa = -^=e ikxx + ^= e - ikix + ^L e ~ ik2X (4) 



<ki y/ki v/c 2 
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FIG. 1: A ring connected with an infinite wire. A (^-potential impurity is at position X. 



Similarly, in region II, III, IV and V we get 
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The lead is along x direction and region III is along y direction, z is used to denote coordinate 
inside the ring, the geometry of the ring being taken care of by applying periodic boundary 
conditions to z. r' u , r' l2 , g'n and g' 12 are the scattering matrix elements and A, B, C, D, E, F, 
G, H , J, K, L, M are to be determined by mode matching. 



Note that at P and Q we have a three legged junction that is schematically shown in Fig. 



2. So far a popularly used scattering matrix for a three-legged junction is 20 
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with a s = |(\A — 2e — 1), 6 S = — 2e + 1) and < e < 0.5. Such junction S'-matrix does 

not include channel mixing and can not account for any contribution from evanescent modes. 



In this work we propose a three-legged junction scattering matrix Sj for a two channel 
quantum wire that can be easily generalized to any number of channels. It is given by 
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FIG. 2: A 3-leg junction. 



Mode matching at the junction P [Fig. [T] gives 
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Similarly, mode matching at the junction Q [Fig. [T] gives 
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Here Zi, Z2 and I3 are shown in Fig. CD a + /5 
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27T0/0O; being the magnetic flux and 



/ic/e is the flux quantum. Eqn. [13] automatically applies periodic boundary conditions to 
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wave function in the ring. Mode matching at the impurity site X [Fig. [T] gives 
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where — . 



— z- 



Tpp' 



(15) 



1 + 2«e * 2fcp 

represents sum over all the evanescent modes and represents sum over all the prop- 
agating modes, p or p' can take values 1 and 2 as there are two propagating modes. 
K e — y ^wr ~ ^?p^> where e = 3,4, ... The inter-mode (i.e. p ^ p') transmission amplitudes are 



t PP ' = f pp i and intra-mode transmission amplitudes are t 
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1 + r pp . Y pp i is given as 



2m7 . r p?r(g + W/2). . r j/7r(g + W/2) . 
— s-sin — sin - 



(16) 



where g is used to denote the position coordinate of the (^-potential impurity. 

We calculate A, B, C, D, E, F, G, H, J, K, L and M numerically from Eqn. 
and Eqn. [TH Persistent current is defined by 



Eqn. 
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which can be simplified to give 
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FIG. 3: I/Iq vs 8ir 2 mEW 2 /h 2 with (a) 7 = (solid line) and (b) 7 = 4 (dashed line). The system 
parameters are l\ = I2 = I3 = 1, a = (3 = 0.3. 

Here 

I {kl) = 2I (\E\ 2 - \F\ 2 + \G\ 2 - \H\ 2 ) {kl) (19) 

is the current when electron is incident along the left lead in k\ channel. This is the scattering 
problem defined by Eqn. [121 Eqn. [13] and Eqn. [TH Similarly, 

I {k2) = 2I (\E\ 2 - \F\ 2 + \G\ 2 - \H\ 2 ) {k2) (20) 

is the current when electron is incident along the left lead in /c 2 channel and this scattering 
problem has to be solved by using a similar set of equations. Here, Iq = 2 mw 2 • 

The nature of current obtained from Eqn. [TH] is shown in Fig [HI We take the energy 
range (4vr 2 < 2m] ™ 2 < 9^2^ i e . 40 < ^ < 88 ) in such a way that both the mo d es are 

propagating. This gives the behavior that is captured in earlier works^*^^. The single ring 
current can be paramagnetic as well as diamagnetic as can be seen from Fig. [3j We shall show 
that when we make one of the modes evanescent, we will get a behavior that is not discussed 
before. As soon as 2mEW 2 /h 2 becomes less than An 2 , the second channel becomes evanescent. 
This is because ki is real, whereas k 2 is imaginary (k 2 —>■ ik 2 in this regime). No electron can 
be incident into an evanescent channel from infinity and so 1^ will not exist. But electrons 
incident in k\ channel can be excited into an evanescent channel in the ring implying that G 
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FIG. 4: I{ l} /I and if 1 ' /Iq vs 8ir 2 mEW 2 /h 2 . The system parameters are h = l 2 = h = 1, 
a = j5 = 0.3, 7 = 4. 

and H in Eqn. [19] are non-zero. A single impurity can excite an electron into the evanescent 
second channel. Scattering at the junctions can also excite an electron into the evanescent 
second channel. An electron residing in an evanescent state will carry a current. This naturally 
arises in the scattering problem that is defined in Eqns. [T21fT4l Evanescent mode current can 
be calculated by directly applying Eqn. [T7] to evanescent mode wave-functions or it can be 
calculated by analytically continuing propagating mode current to below the barrier. Both 
results are equal. 

The S matrix becomes 2x2 and is given by 



Although the S matrix is 2 x 2, its calculation has to be done by using the 6x6 junction 
matrix Sj defined in Eqn. [10] and the 4x4 impurity S-matrix defined in Eqn. [TH g' l2 , r' l2 etc 
are still non-zero, although they do not carry any current but they define the coupling to the 



Unitarity should imply | r' n | 2 + | g f u | 2 = 1 and indeed we get this from the junction matrix 
defined by Sj in Eqn. [TU] and impurity ^-matrix defined in EqnJT4l This implies that Sj is 
appropriate to account for realistic multichannel situations. Sjj does not take into account such 




(21) 



evanescent mode. So Eqns. H2l H3] and H3] still holds with /c 2 — > m 2 where k 2 
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FIG. 5: I/Iq vs a+P when second channel is evanescent. Here l\ = 1.0, h = 0.2, 1% = 1.8 and 
7 = -4.11. 

effects and does not allow one to include coupling to evanescent modes, maintaining unitarity. 
Current is expected to be continuous as the energy changes continuously from evanescent modes 
to propagating modes at 2m ^ /2 = 47r 2 . This also comes out in our calculations 

Note from Eqn. EH that 1^ = if l} + if l} , where if = 2/ (| E \ 2 — \ F | 2 )( fcl \ E and 
F being the wave function amplitudes in the propagating channel and if = 2io(| G | 2 — | 
H | 2 ) (fcl) , G and H being the wave-function amplitudes in the evanescent channel. In Fig. H] 
we have plotted I± and if ^ versus 2mi?r4 72 //i 2 . While if 1 can be positive (diamagnetic) 
as well as negative (paramagnetic), if ^ is seen to be only diamagnetic. if ^ is the current in 
an evanescent channel, and there is a fundamental difference with evanescent channel currents 
in ID. In ID evanescent channel current cannot oscillate with energy because evanescent wave- 
function is not of wave nature^ 2 -. In ID we have to introduce an infinitesimal region of the ring 
where the electron can be propagating (evanescent in the rest of the ring), for the persistent 
current to be oscillating between paramagnetism and diamagnetism 23 . But in the present case 
the second channel is evanescent throughout the length of the ring, its wave-function is not 
of wave nature, and yet can oscillate with Fermi energy. The peaks in if ^ are resonance 
effects due to wave nature of electron wave-function wherein at these energies, the electrons 
can spend a long time in the propagating mode. The impurity also gets a long time to pump 
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more electrons into the evanescent mode. So the evanescent mode current I 2 a l so peaks at 
the same energies where l[ peaks (see Fig. HJ) although the evanescent mode wave-function 
is not of wave nature. The difference between them being that while the peaks in l[ kl ^ can 
be in positive direction (diamagnetic) or in negative direction (paramagnetic), the peaks in 
are, always in the positive direction. As impurity configuration changes, these peaks also 
change randomly. But the peaks in 1% always follow the peaks in l[ . One can also see this 
mathematically. Although the evanescent mode wave-function is not of a wave nature, G and 
H are functions of k% and k 2 , due to the non- locality of quantum mechanics. E and F are also 
functions of k\ and k 2 . While l[ kl ^ will fluctuate around zero value, will fluctuate around a 
certain positive value as disorder configuration changes. Apart from this shift, will follow 
the same rules as l[ kl ^ as far as disorder averaging is concerned. Or more appropriately, I 2 
will follow same averaging rules as conductance that fluctuate with disorder, remaining positive 
all the time. It is much easier to take random values of l 2 and I3 to show this for the average 
current. 

The observable current / = 1^ when second channel is evanescent, is plotted versus flux 
in Fig. The figure shows that when a diamagnetic component is present, the response looks 
like that observed in experiments done by Deblock et aR One can further check the validity 
of our explanation by measuring how the response of the ensemble scales with the number of 
rings in the ensemble. One has to go to a large enough ensemble so that the first harmonic 
has averaged to a flux independent diamagnetic component. This component will scale linearly 
with N, the number of rings present in the ensemble, while the flux dependent part will scale 
as VN. 

Solution of Schrodinger equation in multichannel rings consist of evanescent modes that 
are naturally populated due to scattering. These evanescent modes can carry large persistent 
current that are diamagnetic in nature and are as sensitive to disorder as propagating modes. 
Previous attempts to explain the experimental results on persistent current ignore their con- 
tribution as they were thought to be small and insensitive to disorder. They provide a natural 
explanation for the discrepency between theory and experiments. Future experiments should 
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try to isolate the role of the contributions coming from propagating and evanescent modes. 
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